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Abstract 

It is shown that the Euler hydrodynamics for vortical flows of an ideal 
fluid coincides with the equations of motion of a charged compressible fluid 
moving due to a self-consistent electromagnetic field. Transition to the La- 
grangian description in a new hydrodynamics is equivalent for the original 
Euler equations to the mixed Lagrangian-Eulerian description - the vortex 
line representation (VLR) [0. Due to compressibility of a "new" fluid the 
collapse of vortex lines can happen as the result of breaking (or overturning) 
of vortex lines. It is found that the Navier-Stokes equation in the vortex line 
representation can be reduced to the equation of the diffusive type for the 
Cauchy invariant with the diffusion tensor given by the metric of the VLR. 

PACS: 47.15.Ki, 47.32. Cc 

1. Collapse as a process of a singularity formation in a finite time from the 
initially smooth distribution plays the very important role being considered as one 
of the most effective mechanisms of the energy dissipation. For hydrodynamics of 
incompressible fluids collapse must play also a very essential role. It is well known 
that appearance of singularity in gasodynamics, i.e., in compressible hydrodynam- 
ics, is connected with the phenomenon of breaking that is the physical mechanism 
leading to emergence of shocks. From the point of view of the classical catastrophe 
theory |^ this process is nothing more than the formation of folds. It is completely 
characterized by the mapping corresponding to transition from the Eulerian de- 
scription to the Lagrangian one. Vanishing the Jacobian J of this mapping means 
emergence of a singularity for spatial derivatives of velocity and density of a gas. 
In the incompressible case breaking as intersection of trajectories of Lagrangian 
particles is absent because the Jacobian of the corresponding mapping is fixed, in 
the simplest case equal to unity. By this reason, it would seem that there were 
no any reasons for existence of such phenomenon at all. In spite of this fact, as 
it was shown in @, breaking, however, is possible in this case also. It can 

happen with vortex lines. Unlike the breaking in gasodynamics, the breaking of 
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vortex lines means that one vortex line reaches another vortex line. For smooth 
initial conditions breaking happens first time while touching vortex lines at a single 
point. In the touching point the vorticity becomes infinite. And this is possible in 
spite of incompressibility of both divergence- free fields, i.e., vorticity and velocity. 
To describe the breaking of vortex lines in the papers 0, |^ it was suggested the 
vortex line representation - a mixed Lagrangian-Eulerian description when each 
vortex line is labeled by a two-dimensional marker and another parameter defines 
the vortex line itself. 

This paper is devoted to development of this method to apply to both ideal 
and viscous fiuids. We clarify the role of the Clebsch variables in the vortex line 
representation: these variables can be used as Lagrangian markers of vortex lines. 
However, as well known, these variables can be introduced always only locally and, 
generally speaking, can not be extended for the whole space. In the general situa- 
tion we demonstrate in this paper that transition to the vortex line representation is 
equivalent to consideration of a new compressible hydrodynamics of a charged fiuid 
flowing under action of a self-consistent electromagnetic field. In this case the elec- 
tric and magnetic fields satisfy the Maxwell equations. The most essential property 
of a new hydrodynamics is a compressibility of a new fiuid that for its Lagrangian 
description means compressibility of the corresponding mapping and, respectively, 
a possibility of a breaking. In terms of the Eulerian characteristics this results in 
the breaking of vortex lines when the vorticity ft = curl v takes infinite value. In 
the framework of the new hydrodynamics of a charged fiuid the role of density plays 
the quantity inverse to J which is naturally called as a density of vortex lines. This 
quantity appears from the Cauchy formula for the vorticity f2. Evolution of the 
vortex line density in time and space is defined by the velocity component normal 
to the vorticity. As it is shown in this paper the Cauchy formula can be obtained 
from a "new" Kelvin theorem as well as from the analog of the Weber transforma- 
tion. As the result, the Euler equations turn out to be resolved with respect to the 
Cauchy invariants, i.e., relative to the infinite number of integrals of motion. In 
this case one can consider the Euler equations as the partially integrated equations. 
This circumstance is very important for numerical solution of the Euler equation. 

The vortex line representation can be applied not only to ideal hydrodynamics 
but also to fiow description of viscous incompressible fiuids in the framework of the 
Navier-Stokes equation. In the paper we obtain the equation of the diffusion type 
describing dynamics of the Cauchy invariant in the viscous case with the "diffusion 
tensor" determined by the VLR metric. In its form this equation coincides with 
the equation derived in 0. In this case the equations of motion of vortex lines 
in its original (for ideal fluids) form are understood as the equations given the 
transformation to a new curvilinear system of coordinates. The obtained exact 
equations for description of viscous flows can be considered as the result of exact 
separation of two different temporal scales: the inertial (in fact, nonlinear) scale 
and the viscous one. 
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2. As well known ( see, for instance, 0, []T0[) the Euler equations for an ideal 
incompressible fluid, 

— + (vV)v = -Vp, div V = 0, (1) 

in both two-dimensional and three-dimensional cases possess the infinite (continu- 
ous) number of integrals of motion. These are the so called Cauchy invariants. The 
most simple way to derive the Cauchy invariants is one to use the Kelvin theorem 
about conservation of the velocity circulation. 



/(v-rfl), (2) 



where the integration contour C[r(t)] moves together with a fluid. If in this expres- 
sion one makes a transform from the Eulerian coordinate r to the Lagrangian ones 
a then Eq. can be rewritten as follows: 

oak 

where a new contour C[a] is already immovable. Hence, due to arbitrariness of the 
contour C[a] and using the Stokes formula one can conclude that the quantity 



I = rota ii— (3) 




conserves in time at each point a. This is just the Cauchy invariant. If the La- 
grangian coordinates a in coincide with the initial positions of fluid particles the 
invariant I is equal to the initial vorticity r2o(a). 

Conservation of these invariants, as it was shown first by Salmon [^, is conse- 
quence of the special (infinite) symmetry - the so-called relabeling symmetry. The 
Cauchy invariants characterize the frozenness of the vorticity into fluid. This is a 
very important property according to which fluid (Lagrangian) particles can not 
leave its own vortex line where they were initially. Thus, the Lagrangian particles 
have one independent degree of freedom - motion along vortex line. From another 
side, such a motion as it follows from the equation for the vorticity 

— = rot[vxfi], (4) 

does not change its value. From this point of view a vortex line represents the 
invariant object and therefore it is natural to seek for such a transformation when 
this invariance is seen from the very beginning. Such type of description - the vortex 
line representation - was introduced in the papers 0, 0] by Ruban and the author 
of this paper. 
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3. Consider the vortical flow {ft ^ 0) of an ideal fluid given by the Clebsch 
variables A and /i: 

n = [vxx Vjj]. (5) 

The geometrical meaning of these variables is well known: intersection of two sur- 
faces A = const and /i = const yields the vortex line. It is known also that in 
the incompressible case the Clebsch variables are Lagrangian invariants, being un- 
changed along trajectories of fluid particles: 



dX / , 
- + (vV)A 



0; 



|.(VV), 



0. 



(6) 



Therefore these variables can be taken as markers for vortex lines. It is easily to 
establish that transition in (E) to new variables 



X = X{x,y,z), fi = fi{x,y,z), s = s{x,y,z), 
where s is the parameter given the vortex line, leads to the expression 

1 OR 



where 



is the Jacobian of the mapping 



J 



J ds 

d{x,y,z) 
5(A,/i, s) 



r = R(A, /i, s) 



(7) 

(8) 
(9) 

(10) 



The transform ([T0|) inverse to defines the corresponding transition to the curvi- 
linear, connected with vortex lines, system of coordinates. 

The equations of motion of vortex lines - the equations for R(A, yU, s, t) - can be 
obtained directly from the equation of motion for the vorticity (^). However, the 
most simple way to derive them is to use the combination of the equations (0) : 



V/i 



- VA 



|.(VV). 



0, 



(11) 



which is identical to (|]) due to a linear independence of the vectors VA and V/i. 

Performing in (^) the transformations (J^), we arrive at the equation of motion 
for vortex lines 0: 



dR 

ds 



dR 

'dt 



v(R,t) 



0. 



(12) 



This equation has one important property: any motion along a vortex line does not 
change the line itself. It is easily to check that Eq. (|12]) is equivalent to the equation 



— = v„,(R,t), 



(13) 
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where v„ is the velocity component normal to the vorticity vector. 

In accordance with the Darboux theorem, the Clebsch variables can be intro- 
duced locally always but not globally. It is well known also that the flows pa- 
rameterized by the Clebsch variables has a zero helicity integral / (v ■ rot v)(ir - 
the topological invariant which characterizes a degree of knottiness of vortex lines. 
Therefore to introduce the vortex line representation for flows with nontrivial topol- 
ogy it is necessary to come back to the original equations of motion (|l|) and @ for 
velocity and vorticity. 

4. According to the equation the tangent to the vector f2 velocity component 
does not effect (directly) on the vorticity dynamics, i.e., in (|D we can put, instead 

of V, its transverse component v„. 

The equation of motion for the transverse velocity v„ follows directly from the 

equation (|l]). It has the form of the equation of motion of charged particle moving 

in an electromagnetic field: 

-g^ + (v„V)v„ = E + [v„ X H], (14) 

where the effective electric and magnetic fields are given by the expressions: 

^^r, 




dt ' 



(15) 



H = rot v^. (16) 

Interesting to note that the electric and magnetic fields introduced above are ex- 
pressed through the scalar ip and vector A potentials by the standard way: 

^=p+^, A = v,, (17) 
so that two Maxwell equations 

dH 

div H = 0, — = -rot E 
dt 

satisfy automatically. In this case the vector potential A has the gauge 

div A = —div v„, 

which is equivalent to the condition div v = 0. 

Two other Maxwell equations can be written also but they can be considered as 
definition of the charge density p and the current j which follow from the relations 
(|I5D and (^6|). The basic equation in the new hydrodynamics is the equation of 
motion (|14|) for the normal component of the velocity which represents the equation 
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of motion for nonrelativistic particle with a charge and a mass equal to unity, the 
light velocity in this units is equal to 1. 

The equation of motion (|T^) is written in the Eulerian representation. To trans- 
fer to its Lagrangian formulation one needs to consider the equations for "trajecto- 
ries" given by the velocity v^: 

^ = v„(R,t) (18) 

with initial conditions 

R|t=o = a. 

Solution of the equation ([T8| ) yields the mapping 

r = R(a,t), (19) 

which defines transition from the Eulerian description to a new Lagrangian one. 
The equations of motion in new variables are the Hamilton equations: 

^ dh ■ dh , , 

R-=^, (20) 



an' dP' 

where dot means differentiation with respect to time for fixed a, P = v„ + A = v 
is the generalized momentum, and the Hamiltonian of a particle h being a function 
of momentum P and coordinate R is given by the standard expression: 

i.e., coincides with the Bernoulli "invariant". 

The first equation of the system (|20| ) is the equation of motion ([1^), written in 
terms of a and t, and the second equation coincides with ([TB|). 

For new hydrodynamics (0) or for its Hamilton version (pO]) it is possible to 
formulate a "new" Kelvin theorem (it is also the Liouville theorem): 



r = ^(p-rfR), (21) 



where integration is taken along a loop moving together with the "fluid". Hence, 
analogously as it was made before while derivation of (|^) we get the expression for 
a new Cauchy invariant: 

I^.ot„(p,g. (22) 

Its difference from the original Cauchy invariant (|^) consists in that in the equation 
of motion (p!8| ) instead of the velocity v stands its normal component v„. As conse- 
quence, the "new" hydrodynamics becomes compressible: div v„ ^ 0. Therefore on 
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the Jacobian J of the mapping (^) there are imposed no restrictions. The Jacobian 
J can take arbitrary values. 

From the formula (p2D it is easily to get the expression for the vorticity f2 in the 
given point r at the instant t (compare with 0, 0]): 

n(..i) = <"W-^;)^<°'". (23) 

where J is the Jacobian of the mapping (^) equal to 

9(Xi,X2,X3) 



J 



(9(ai,a2,a3) 



Here we took into account that the generalized momentum P coincides with the 
velocity v, including the moment of time t = 0: Po(a) = vo(a). r2o(a) in this 
relation is the "new" Cauchy invariant with zero divergence: djva^oio) = 0. 

The representation (|23|) generalizes the relation to an arbitrary topology of 
vortex lines. The variables a in this expression can be considered locally as a set of 
A, n and s. 

As known (see, for instance, 0]), expression for the Cauchy invariant can be 
obtained from the Weber transformation. This is the representation of velocity in 
terms of the initial data which can be obtained by integrating the Cauchy invariant 

Consider the following one-form uj = {P ■ dR) and calculate its time derivative. 
By means of the equations of motion (|20D we get: 



Cj = d[-h + (PR)]. 
Hence it follows that the vector function 

Uk = ^ Pi, 

dak 

dependent on t and a, will obey the following equation of motion: 

Integration of this equation in time gives the Weber-type transformation: 

Mfc(a, t) = Mfco(a) + (24) 
oak 

where the potential $ satisfies the nonstationary Bernoulli equation: 

^22 
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If $|t=o = then the time independent vector uo(a) coincides with the initial 
velocity vo(a). By applying the operator curl to the relation (j2^) we arrive again 
at the Cauchy invariant p^). 



Thus, in the general situation the equation of motion of vortex lines has the 
form (13) which is completed by the relation (p3|) and the equation 



0(r,t) = rotrv(r,t) (25) 

with additional constraint divrv(r,t) = 0. 

The equations of motion ([l8|) , (p5|) together with the relation ( p3l) can be con- 
sidered as the result of partial integration of the Euler equation (|l]). These new 
equations are resolved with respect to the Cauchy invariants - an infinite number 
of integrals of motion, that is a very important issue for numerical integration (see 
[D, H). For the partially integrated system the Cauchy invariants conserve auto- 
matically that, however, for direct numerical integration of the Euler equation one 
needs to test in which extent these invariants remain constant. Probably, this is one 
of the main restrictions defining accuracy of discrete algorithms of direct integration 
of the Euler equations. 

Another very important property of the vortex line representation is absence 
of any restrictions on the value of the Jacobian J which do exist, for instance, for 
transition from the Eulerian description to the Lagrangian one in the original Euler 
equation (|1|) when Jacobian in the simplest situation is equal to unity. The value 
1/J for the system (|18|), (^51) , (p3D has a meaning of a density n of vortex lines. 



This quantity as a function of r and t , according to (p!8| , obeys the discontinuity 
equation: 

— -fdiv,(nv„) = 0. (26) 

In this equation div^v^ ^ because only the total velocity has zero divergence. 

5. Consider now the question about application of the VLR to flows of viscous 
fluids. Write down the Navier-Stokes equation for vorticity fi: 

— = rot[v X r2] - z/ rot rot fi, (27) 

and perform in this equation the transform to new variables a and t by means of 
changes defined by the equation ([TsD together with the Cauchy relation (^3]) where 
Oq is assumed a function of not only a but also time t: Oq = fio(ci, t)- 

Then after substitution (^31) into (|27|) the first term in the right hand side is 
cancelled because of (|1^). At the result, the equation ( p?] ) is written in the form: 



J \~dr 



vJR=-^^rotrot|^(^^o■V,)R|. (28) 



Next, change differentiation relative to r in the r.h.s. of (p8[) to differentiation 
against a. After simple, but cumbersome calculations the equation (123) transforms 



8 



into the equation for fio(a, t): 

= )• 

Formally it is a linear equation for fig, here g is the metric tensor equal to 

dxi dxi 

gal3 = ^ ■ ^ • 

oaa oap 

The equation ( p9|) for the Cauchy invariant formaly coincides with that ob- 
tained by Zenkovich and Yakubovich for incompressible hydrodynamics |P in which 
the variables a are assumed to be Lagrangian markers of fluid particles. In the 
Zenkovich- Yakubovich equation the Jacobian J is proposed to be independent on 
time, in the simplest case equal to 1. Just this is a principle difference between the 
Zenkovich- Yakubovich equation and the equation (^). J in ( p9|) is a function of 
time t and coordinates a. 

Remarkable peculiarity of the obtained system is the exact separation of two dif- 
ferent temporal scales, responsible for the inertial (in fact, nonlinear) processes and 
for the viscous processes. The former ones are described by the equation (0), and 
the latter by the equation of diffusive type (^) in which the diffusion " coefficient" , 
proportional to viscosity z/, is defined by the metric of the mapping r = R(a, t). 
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the paper [§] before its publication. This work was supported by the RFBR (grant 
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